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ABSTRACT 

In 1969 Edward Conklin measured the anisotropy in celestial emission at 8 GHz with a resolution 
of 16.2° and used the data to report a detection of the CMB dipole. Given the paucity of 8 GHz 
observations over large angular scales and the clear evidence for non-power law Galactic emission near 
8 GHz, a new analysis of Conklin's data is informative. In this paper we compare Conklin's data 
to that from Haslam et al. (0.4 GHz), Reich and Reich (1.4 GHz), and WMAP (23-94 GHz). We 
show that the spectral index between Conklin's data and the 23 GHz WMAP data is /3 = —1.7 ± 0.1, 
where we model the emission temperature as T oc v$ . Free-free emission has j3 rj —2.15, synchrotron 
emission has /3 ~ —2.7 to —3. Thermal dust emission (f3 « 1.7) is negligible at 8 GHz. We conclude 
that there must be another distinct non-power law component of diffuse foreground emission that 
emits near 10 GHz, consistent with other observations in this frequency range. By comparing to the 
full complement of data sets, we show that a model with an anomalous emission component, assumed 
to be spinning dust, is preferred over a model without spinning dust at 5tx (A% 2 = 31). However, the 
source of the new component cannot be determined uniquely. 



1. INTRODUCTION 



Many observations have shown that there is diffuse 
dust-correlated emission at frequencies below 40 GHz 
where the thermal emission by dust grains is usually ex- 
pected to be negligible. The effect was first seen at large 
angula r scales in the COBE DMR data by |Kogut et al. 



co-location of thermal dust and free-free emission. These 
observations were confirmed b y ground based measure 



ments at finer angular scales (de Oliveira-Costa et al 



1997| |Leitch et al.|1997[) . Soon thereafter, it was realized 



that the dust-correlated emission might in fa ct be due to 
microwave emission by sp inning dust grains ( Jones|1997 
Draine & Lazarian 19981). In other words, where there 



is more thermally emitting dust radiating at frequencies 
v > 100 GHz, there is also more spinning dust emit- 
ting near v ~ 20 GHz. Since then, the correlated diffuse 
emissi on has been seen at high and low Galactic lati- 
tudes ([de Oliveira-Costa et al.||1998l |1999| |Leitch et al. 
2000} [MyTKerjee et al.||^001MHamilton U Gangal|Ml| 
Mukherjee et al ||2002[ |2U05f |Beniiett et al.|J2003[ |Ban- 
day et al.j|2003| |Lagache | |2003| |de Oliveira-Costa et al 



2004| |Finkbeiner||2004[ |Davies et~al 
Cerezo et al 



et al 
2009 



2006[ |Fe rnandez- 
20061 |Boughn fc Pober|[20071 Hildebrandt 
20071 IDobler fc Finkbeiner||2008a| b[ I D obler et al. 

sard 

Kogut et al.l I2011| [Gold et ;il.||201 \\. At large' angu- 
lar scales, WMAF ( jbennett et al.||2003[ ) suggested that 
this correlation could also be explained by variable- index 
synchrotron emission co-located with dusty star-forming 
regions. However, alternative and subsequent analyses 
gave more support to the spinning dust hypothesis and 
attributed the source of t he correlation to a population 
of small spinning grains (|de Oli veira-Costa et a l. | [2002j 
Lagache||2003[ |Ysard et al.||2010[ |Macellari et al.||201ip . 
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(1996), who attributed the phenomenon to the simple fomn^rffT^ 



Dust-correlated emission at microwave frequencies has 
also been observed in more compact region s (e.g ., 
i Finkbeiner et all (f2002l |2004[>; [Watson et all (|2005j); 
Casassus ct al. (2004); Battistclli et al. (2006); Casas- 
sus et a l.| (|2006|); |lglesias-(jroth| ( |200 6); Casassus et al. 
Dickinson et al.| J2009D; |Scaite et al.l (|2009[); 

1'ibbs et all (|2010[); jScaife et at 
lckinson et al.l (2010); Planck Collab oration 
as! cllanus ci al.||2 



(2008) 



iv et al 
2010a b); DT 
et al.|ti20llt 



ms on et al 



Vidal et al (201 
Lopez-Caraballo et al.| ( |2011[ ) 



(2010) 

iTT mi 



Most recently, the 



2011) 
'lanck 



satellite has shown in detail the presence of a non power- 
law component. A new component of emissi on is clearly 



seen in the Perseus and p Ophiuchus regions ( |Planck Col 
laboration et al.||2011 >. The component is well modeled 



as spinning dust, although the range of gas temperatures 
is large. 

It has long been realized that large area maps in the 
5-20 GHz range would be ideal for separating free-free 
from variable-index synchrotron emission, spinning dust, 
or any other emission process. For diff use sources, COS- 
MOSOMAS tiHildebrandt et al.||2007"| observed the sky 
between 11 and 17 GHz and ide ntified a dust-correlated 
comp onent, as did TENERIFE ( |de Oliveira-Costa et al. 
1999| ) obser ving at 10 and 15 GHz. The ARCADE l l 
experiment ( Kogut et al.||201l" > observed at 3, 8, and 10 



GHz and also found evidence m support of a ne w compo 
nent. Here, we show that Conklin's 1969 data (Conklin 



1969a|b[ ) further improve constraints in this frequency 



range and at large angular scales. 

Although spinning dust is currently the preferred sin- 
gle hypothesis for explaining the correlation, there are 
10 free parameters in the spinning dust model (e.g. 



Ali-Haimoud et al. (20091); and there could b e multiple 



processes, inclu ding magnetized dust emission ( Draine & 
Lazarian][L99 9 ) , at work 



The rest ot the paper is outlined as follows. We de- 
scribe the observations in <f2l discuss an interpretation of 
the observations in fj3j andconclude in Sj4j 
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2. CONKLIN'S OBSERVATIONS 

Conklin's observations took place on White Mountain 
in 1968 and 1969, using a coherent receiver at 8 GHz 
with an effective full width at half maximum beam of 
$1/2 — 16.2°. The radiometer Dicke-switched at 37 Hz 
between two 9^ 2 — 14.5° feeds pointed ±30° from the 

zenith along an E-W baseline. The data were averaged 
over 4 minutes, the entire apparatus was rotated 180° 
over one minute, data were averaged again for 4 more 
minutes, and finally the apparatus was rotated back to 
its original position. During each such ten-minute cy- 
cle, the temperature difference between the two feeds 
was recorded; Conklin reported the differences as "east"- 
"west." The 10-min differences were then averaged with 
a 30-min FWHM Gaussian function, effectively broaden- 
ing the beam profile to 16.2° along the scan direction. 
The data we use, shown in Figure [I] and reported in 
Table nl come from the Nature publication and Con- 
klin's Ph.D. thesis. For this analysis, we use Conklin's 
raw data after subtracting the contribution from the now 
well-established dipole. The per-point error bars are de- 
duced from the "probable error" reported in his thesis. 
They may be treated as independent. They are larger 
than the purely statistical error by a factor of 1.6 R 

The measurements were made at a single celestial cir- 
cle at S — 32° as indicated in Figure [2] The highest 
amplitude point at RA= 300° comes from when the "E" 
beam crosses the Galactic plane at b = 0° and I = 69.5°, 
and the "W" beam is out of the plane. The "E" beam 
again crosses the Galactic plane at I = 176.4°. 

After subtracting up to a 10 mK Galactic contri- 
bution by extrap olating the 404 MHz Pauliny-Toth & 
Shakeshaft map (Pauliny-Toth & Shakeshaft 19621 to 
8 GHz, Conklin reported a preliminary measurement of 
the dipole with amplitude 1.6 ± 0.8 mK in the direction 
a = 13h. After Conklin's 1969 Ph.D. thesis, this be- 
came 2.3 ± 0.9 mK in the d irection a = l lh as reported 
at the IAU 44 symposium ( |Conklin||l9"72"l ). Though the 
error bars include an estimate tor a ±0.1 uncertainty 
in the Galactic index (the statistical uncertainty was 
« 0.02 mK) , there were li ngering doubts about the ex- 
trapolation ( Webster|1974 ) . Neglecting the correction for 
the Earth's motion, the WMAP dipole (3.358±0.017 mK 
in direction a = 11.19h at 8 — —6? 90 for the full sky) has 
amplitude 2.83 mK at 8 = 32°, in excellent agreement. 

2.1. Calibration and beam characteristics 

The system was calibrated with a IK argon noise tube 
that injected power into the feeds through a cross guide 
coupler. This is effectively a full-beam calibrator. Con- 
klin was mindful of the required stability and anchored 
all components with thermal straps every two inches. He 
measured the losses and reflections from all components 
and accounted for the atmosphere. A formal error is not 
reported but the system is similar t o the one he used for a 
10.7 G Hz anisotropy measurement ( Conklin & Bracewell 
1967) where he reported ±10%, which we adopt. 

Understanding the beam profile is critical for compar- 

3 To ensure accuracy in calculations, we digitized Conklin's data 
at twice the resolution reported in Table [TJ using 48 bins in RA 
with error bars increased by \/2. For the results we report, the 
calibration uncertainty dominates the errors in the Table. 
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Fig. 1. — Conklin's observations of the anisotropy at 8 GHz. 
The data are plotted for each differential E-W observation at S = 
32°, with the now well-established dipole subtracted. The dashed 
line is the WMAP dipole. The offset is not constrained by the 
measurement. The peak at RA ~ 300° corresponds to the "E" 
beam crossing the Galactic plane at b = 0° and I = 69.5°. 



TAB LE 1 

Data at 8GHz from Conklin ( 1969a|b 



RA 
(deg) 


Total 
anisotropy (mK) a 


Dipole- 
removed (mK) b 





24.4 ± 1.13 


4.38 


15 


19.8 ± 1.13 


0.60 


30 


16.5 ±1.13 


-1.94 


45 


15.1 ± 1.13 


-2.86 


60 


16.6 ±1.13 


-1.00 


75 


18.8 ± 1.13 


1.46 


90 


21.3 ± 1.13 


3.81 


105 


23.4 ± 1.13 


5.65 


120 


24.7 ±1.27 


6.42 


135 


23.1 ± 1.27 


4.20 


150 


22.1 ± 1.27 


2.39 


165 


21.7 ± 1.27 


1.13 


180 


22.1 ± 1.27 


0.65 


195 


22.7 ± 1.27 


0.50 


210 


23.0 ± 1.27 


0.10 


225 


22.2 ± 1.27 


-1.31 


240 


19.4 ± 1.27 


-4.50 


255 


12.4 ± 1.27 


-11.6 


270 


5.67 ± 1.27 


-18.3 


285 


14.5 ± 1.27 


-9.13 


300 


20.4 ± 1.13 


-2.82 


315 


23.8 ± 1.13 


1.30 


330 


29.3 ± 1.13 


7.58 


345 


34.2 ± 1.13 


13.3 



a The values include an offset that is removed in our analysis and 
in Figure 1. 

b The errors are the same as for the total anisotropy. Here the 
offset has been subtracted. 



ing to other data to assess the Galactic contribution. 
During observations, the sky drifts through the H-plane, 
the profile of which was measured. A measurement of 
the E-plane was not done. Based on his dimensions, 
we recomputed the full beam pattern using a more pre- 
cise calculation ( Sletten||l988 ). We find an H-plane pat- 
tern 6^ 2 = 14.5°, in agreement with Conklin, and an 



E-plane pattern with 9f^ 2 = 13.7° as compared to Con- 
klin's 9^ 2 = 14°. For the full beam, we compute a 6% 



8 GHz emission 



3 
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Fig. 2. — Conklin's observations with respect to the Galaxy. The 
"E" and "W" beams are separated by 60° and are shown here as 
circles of diameter 16.2°. The dashed line marks the Galactic plane 
and the grey swath indicates b = ±15°. 

smaller solid angle. This translates into a possible bias 
that the peak-to-peak amplitude of the data is low by 
6%. We account for this possibility in the analysis. 

3. COMPARISON TO OTHER DATA SETS 

We compare Conklin's data to that from Haslam 
dHa slam et al.||19811 11982[) at 0.408 GHz; Reich and Re 
Reich & Reich|1986 



ich (Reich 
at OTTJ 



1982; 



mrWMAP at 23, 33. 



Reich et al. 



2001 
and 94 G 



'If 

)H; 



( Jarosik et al. 2011 1; and the FDS map mo del 8 extrap- 
olated to 94 GHz pinkbeiner et al.||1999[) . We adopt 
calibration errors ot 7%, 4%, and 2% tor Haslam, Reich 
& Reich, and WMAP respectively. The statistical errors 
are negligible compared to the uncertainty on Conklin's 
observations. We convolve the multi-frequency maps to 
the same 16.2° resolution and extract differenced mea- 
surements at 5 = 32° to compare to the Conklin data, 
shown in Figure HI We do not account for the asymmetry 
in the Conklin effective beam, but test that varying the 
beam size by 2° degrees has a negligible effect on con- 
clusions. In all of our fits the calibration uncertainties 
dominate the statistical uncertainties. 

We initially compare the emission, /, at Vq = 8 GHz 
to the WMAP K-band emission (v = 23 GHz), using a 
two parameter model, 



I(yo,n) = I{v,n){v /vf + A. 



(1) 



We fit for the index, (3, and offset, A, finding a marginal- 
ized limit of £ = -1.7 ± 0.1, with best-fit x 2 = 64 (re- 
duced x 2 = 1-34). Because WMAP and Conklin are 
not observing the same mixture of foreground compo- 
nents, the high x 2 is n °t surprising. Figure [3] shows 
the marginalized distribution for the spectral index of 
the model. We include calibration error and the pos- 
sible 6% calibration bias. The index is incompatible 
with free- free emission (/3 = —2.15), synchrotron emis- 
sion —3 < P < —2.7, or a "breaking" synchrotron index. 
Additionally, the WMAP best-fitting maximum entropy 
(MEM) model, excluding a spinning dust component, 




-2.0 

Spectral index 

Fig. 3. — Marginalized probability distribution for the effective 
spectral index between the Conklin 8 GHz data and the WMAP 
23 GHz K-band data. The best fit model has \ 2 = 64.5 for 46 
degrees of freedom. The index is shallower than would be expected 
from sychrotron (left vertical dashed line) and free-free emission 
(right vertical dashed line), indicating the need for an additional 
component. The curves are, from left to right, Conklin plus a 
6% systematic bias and calibration error, nominal Conklin with 
calibration error, and nominal Conklin without calibration error. 
The leftmost index shown is f) = —1.72 ± 0.1. 



TABLE 2 

Scaling factor between Conklin's data and other 
data sets. 



Data Frcq (GHz) 


s(u) x 2 /DOF 


Haslam 0.408 
Reich 1.42 
Conklin 8 
WMAP 23 
WMAP 33 
WMAP 41 
WMAP 61 
WMAP 94 


2890 ± 130 2.04 
107 ±4.54 1.75 

1 

0.179 ±0.0075 1.33 
0.0732 ± 0.0031 1.54 
0.0431 ± 0.0018 1.81 
0.0194 ± 0.00086 2.69 
0.0158±0.00072 3.19 



4 Obtained from de O liveira-Costa's Global Sky Model 
|Qliveira-Costa et al.|2008[ l 



:de 



over-predicts the emission at 8 GHz by a factor of two. 
We conclude that there must exist another component of 
foreground emission. 

To determine the scaling between Conklin's data and 
each of the data sets in the range ^=0.4 — 94 GHz, we fit 
I{v) = s(y)I(yo) + A(y). Figure 111 shows how the inten- 
sity measured at each frequency, l{v), compared to the 
scaled Conklin data, s{v)I(vq) + A{v\. The estimated 
scaling factors are reported in Table [2] with the good- 
ness of fit. The closest fit is obtainea between 8 GHz 
and the WMAP K band because these are the closest 
in the logarithm of the observation frequency. A sim- 
ple power law extrapolation is a poor fit for more widely 
spaced frequencies as multiple components contribute to 
the emission. Figure [5] shows the frequency dependence 
of the estimated scaling factors s(v) after multiplying by 
the rms of Conklin's data. Consistent results are found 
by fitting to the I = 176° or I — 69.5° Galaxy crossings, 
the rms of each data set, or to the peak-to-peak ampli- 
tude as a function of frequency. The fit is dominated by 
the Galactic crossing regions. 
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Fig. 4. — Emission intensity at frequencies 0.4 < v < 94 GHz 
in the 8 = 32° strip surveyed by Conklin compared to the scaled 
uq = 8 GHz signal. The data from each map, I(y), is fit to the 
Conklin data using I(u) = s(v)I(vq) + A(y), and compared to the 
scaled Conklin data using the best-fitting s{v) and A(v). With 
multiple emission components, a global scaling of the 8 GHz data 
does not fully describe the foreground intensity. 



To find a physical model that fits Conklin's observa- 
tions as well as the Haslam, Reich and Reich, and WMAP 
data we perform a simultaneous fit to s(v) in the range 
0.4 — 94 GHz. The model includes free-free, synchrotron, 
thermal dust, and spinning dust. We test the improve- 
ment in x 2 as we add a spinning dust component, using 



the model of Ali-Haimoud et al. (2009). As free parame- 
ters, we take the total hydrogen number density riH and 
the gas temperature T, assuming a single density and 
temperature can be used to model the global emission. 
For the other parameter s, we use the values given by 
Ali-Haimoud et al. ( 2009 ) in their sample spectrum and 
nominal parameters m Weingartner & Draine ( 2001 1 Ta- 
ble 1, line 7. Best fits are obtained with «h = 20 cm -3 
and T = 300K. The values of these parameters, however, 
would b e affected by m ore refined models of the emission 
process ( |Draine||2011 ). For the thermal dust emission, 
we consider two models: one in which the thermal dust 
amplitude is constrained by the FDS99 model, and one 
in which the amplitude is allowed to vary (hereafter re- 
ferred to as "free dust" ) . The results are not significantly 
different, so we report only the "free dust" results)^] 

Model 1. First, we fit s{v) with synchrotron, free-free, 
and thermal dust emission: 



Ah{v) = a v ai + a 2 v Z - Lj + a 3 u 



.1-7 



(2) 



where aov ai represents synchrotron emission with spec- 
tral index ai and amplitude eto, <22^~ 215 represents free- 
free emission, and 03 v 7 represents thermal dust emis- 
sion with fixed spectral index. The fit parameters do, a 2 , 
and 0,3 are required to be positive, and the synchrotron 
spectral index a\ is varied in the range —3.5 < a% < —2. 
The probability distribution for the four parameters is 
estimated using Markov Chain Monte Carlo sampling 
methods. The best-fitting synchrotron index and \ 2 are 
reported in Table [3] the model is a poor fit to the data 
with x 2 /dof = 49.9/4. 

Model 2. We then include an additional spinning dust 
component: 



M 2 {v) = a n v ai + a 2 i>- z - L0 + a 3 f L ' + a 4 £> s (f) (3) 

where D s (v) represents the spinning dust emission tem- 
plate, with amplitude 04. Adding the spinning dust im- 
proves the goodness of fit by A% 2 = 31 indicating a 
strong preference for this additional component. 

This model is physically reasonable for a single pixel 
or direction on the sky. Generally, the relative ampli- 
tudes of each emission mechanism spatially varies at any 
given frequency, and the frequency dependence of the 
synchrotron, spinning dust, and thermal dust also has 
spatial variation. In our case, the model assumes that 
the relative foreground intensity over the region of sky 
observed by Conklin may be quantified by a single scaling 
factor s(v). This is likely over-simplified, but the fit is 
good. It is noteworthy that the amplitude of the best fit 
synchrotron spectrum at WMAP's 23 GHz band is low 
compared to the free- free and spinning dust, contributing 
only ~ 10% of the anisotropy, although the best fit syn- 
chrotron spectral index, f3 w —3.00, is reasonable. One 
must keep in mind that with so little data the parame- 
ter degeneracies are large, and a more realistic model of 



5 For all the models described hereafter, 
so that all values are unitless. 
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Fig. 5. — Emission as a function of frequency in the 5 = 32° 
strip, for observations in the range 0.4 < v < 94 GHz relative to 
the 8 GHz observations. The emission is modeled as synchrotron, 
free-free, thermal dust, and spinning dust, with x 2 /dof = 12.8/3, 
compared to 49.9/4 when excluding spinning dust. 

the relative contribution of the foreground components 
would all ow for spatially varyin g amplitude ratios. For 



example, Macellari et al. (2011) estimate ~ 1/3 of the 



all-sky intensity anisotropy to be synchrotron at 23 GHz 

4. DISCUSSION 

Anisotropy in the diffuse microwave emission at 8 GHz, 
originally measured by Conklin in 1969 to estimate the 
CMB dipole, now sheds light on Galactic emission. The 



observed signal is consistent with synchrotron, free-free, 
and spinning dust emission from the Galaxy. In combina- 
tion with multi-frequency observations, the 8 GHz data 
strongly disfavor an emission model with no anomalous 
dust component, assumed to be rapidly rotating PAH 
grains. The presence of this additional diffuse component 
is consistent with observations by ARCADE-2, COSMO- 
SOMAS, and the TENEFJFE experiment, and with tar- 
geted measurements of dusty regions in the Galaxy. Ac- 
curate characterization of the Galactic foregrounds in in- 
tensity and polarization is important for extracting cos- 
mological information from the CMB, and will be vital 
for constraining inflationary models via the large-scale 
polarization signal. Upcoming low-frequency observa- 
tions, for example from the C-BASS experiment at 5 GHz 



(King et al. 



2010), will shed further light on the diffuse 



anomalous dust behavior and allow its properties to be 
better established. 
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